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We considered the phenomenon of the Andreev reflection of quarks at the interface between the 
cold quark-gluon plasma phase and 2SC color superconductor. 



I. INTRODUCTION 

CN ■ 

. The very interesting phenomenon arises at the junction between a normal-metal and a superconductor. The 
' reflection of a conduction electron from the superconducting barrier, first elaborated in the paper is known in the 
literature as the Andreev reflection, ft was shown that the excitation in the normal metal outgoing from the reflection 
^ ■ point has the sign of the momentum, the sign of the charge, and the sign of the effective mass opposite to that of 
^_ ' the incoming electron (hole quantum numbers). This peculiar prediction was checked by tunneling experiments and 



(N 



ballistic (direct) observations 

In this letter we raise the interesting question of the interaction between the phase containing the free quarks (QGP) 
and the color-superconductor (CS). We shall show that the free quarks from the QGP phase falling at the CS barrier 
fS| ' are reflected in the similar fashion as was described by Andreev in 

J> \ In the next section we develop the Bogolubov - de Gennes equations for superconducting phase. The third section 
contains the detailed discussion of the Andreev reflection. In the last chapter we make some comments about the 
^ influence of the Andreev reflection on the phenomena that can take place inside the Neutron Stars. 

^ \ II. BOGOLUBOV - DE GENNES EQUATIONS FOR 2SC PHASE 

O ; 

(-H , Our interest is the description of the phenomena that happens at the QGP/CS junction at moderate densities 
i' expected in the Neutron Stars. In this work we do not address the question of the gluon interaction with the CS 
' , phase thus we just concentrate on the fermionic degrees of freedom. The model of QGP is the gas of the free quarks. 

The quarks in QGP built the Fermi ladder, with the Fermi energy /i (of order 400 — 500 MeV in the Neutron Stars), 
(-H I where is a quark chemical potential. In the case of the two light quarks, which we consider here, the CS phase 
• creates the 2SC superconductor (for the review of the subject see |3j). The essential physics of 2SC phase at moderate 
.5^ I densities, is described by the effective four-fermion point interaction with the attractive pseudoscalar channel Q. This 
. model leads to the Cooper instability of the Fermi sea which creates the new vacuum of condensed Cooper pairs. 
$^ ' The physical excitations are fermionic quasiparticles separated from the vacuum by the energy gap |A|. The energy 
5^ , gap was calculated in the mean-field approximation in the variety of models giving the values 40 — 140 MeV which is 
about 10-30 per cent of the Fermi energy. 

The effective hamiltonian at the mean-field level takes the form: 



H ^ I (fx 



. V + m7o - A*) + y^it(r2),feio/3C75^^^ - ^i^i{r2)Mo.pC^5^0 + ^(A, m)| (1) 

where m is a quark mass and ip^ is the field operator of the quark. The greek indices describe the color and the latin 
indices describe the isospin quantum numbers. The matrix T2 is the antisymmetric Pauli matrix and t is one of the 
three antisymmetric Gell-Mann matrices. For our calculations we chose it arbitrarily as the A2. The gap parameter 
A(f) depends on the position. It vanishes inside the QGP phase and it takes a nonzero value in the 2SC phase. The 
equation of motion describing our system takes the form: 
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z^^ = (-za . V + m7o - m) + A{r){T2)itJ,Cj5^'^^^ (2) 

The equations ^ describe the interaction of the 6 "famiUes" of quarks coupled through the non-zero value of the gap 
parameter. Using standard conventions (e.g. [^), one can find that the up, "red" quark couples to the down "green" 
quark, the up "green" to the down "red", whereas the "blue" quarks remain free. For the first "family" the equations 
(|) read: 

z^R = ■V + mjo-f^)rR- A(f)C75V'G^^ (3) 

ii^'J = -V^^a - 4^ (m7o - /i) - A(f)>]^^C75 

where u and d are isospin and R and G color degrees of freedom. The equations for the second "family" (up "green" 
and down "red" quarks) can be evaluate exactly the same way and the "blue" quarks are not interesting in the first 
approximation. One can change the operator equations to c-number equations by taking the expectation values 
/ = (^'olV'i^l^i) and = where |\l/o) , l^l/i) are ground and excited state respectively. In the QGP 

phase (|A| = 0) the fields / and describe a "particle" and a "hole". Let us now find the quasiparticle wavefunctions 
in the case of A = const. It is convenient to decompose the fermi fields in the plane wave basis. Remembering 
that the relativistic effects are suppressed in the presence of the Fermi sea one can neglect the contribution from the 
antiquarks: 

f{t,r) = ^Usiq)as{q)exp{-iEt + iq-r} (4) 

S 

9\t, r) = J2 dli-'l)K{Q) eM-^Et + iq- r), 

s 

where s describes the spin degrees of freedom, as{q) and Psiq) are some c- numbers. The energy E describes the 
energy of the quark with respect to the Fermi surface in the absence of the gap parameter A. The Dirac bispinors Us 
and ds satisfy the equations: 

{a-q + m7o - ^J.)us (q) = egUs {q} (5) 
dl{-q){a ■ q- m7o + m) = dl{-Q)^q, 

and because momentum q can be in general the complex number one has to carefully distinguish between the left and 
the right eigenvectors. Inserting (^) into the equations (H), using bispinors algebraic relations and assuming constant 
value of the gap parameter one arrives at the Bogolubov - de Gennes equations for a and (3 parameters: 

Ea^{q) = ega^{q) + APl{q), (6) 
E(3liq) = -egPliq) + A*a^{q) 

and the similar set of equations with reverse spin projection on the momentum q: 

Eai{q) = e,aiiq)-Ap^q), (7) 
E(3*^{q} = -egl3*^iq)-A*ai{q) 

where — y/ q^ + — ^ = —iq. We have two possible solutions of the uniform equations (||) (the case of (^ is 

similar): 

for = ^ 




and for 
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"T = Wo exp i- (9) 
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where ^ = ^/ E"^ — | Ap and (5 is a phase of the gap. Let us point out that only for > |A| the solution describes 
the propagating excitations. Finally the wavefunction describing the quasiparticle excitation of energy E in the 2SC 
phase takes the form: 



exp {—lEt + iq2 ■ r) 



exp{-tl)^/^(cdl^{-q-2)-Dd\^{-q2: 



where qi 2 = ±-y/ ± — and A, B, C and D are arbitrary constants. 

III. ANDREEV REFLECTION 

Let us consider the simple physical problem of the quark scattering on the superconducting plane surface placed at 
z=0 in space. If we approximate the boundary as the step function (A = for z < and A = const for z > 0) then 
we are looking for stationary solution of equations (|^), independently for negative and positive z, supplemented with 
the condition that matches the wavefunction at z = 0. 

Let us consider in detail the situation when the "red" quark, with energy E and spin projection f falls at the plane 
boundary from the left. Then the wavefunctions (for A = 0) take the form: 

/ lit (fc) exp ifcz + i^ut (— fc) exp — i/cz \ , , , , 

^<(^'^)-( Hd\-i-,)%L )exp(-^M) (11) 

The coefRcients F and H describes the amplitude of the reflection of the w-particle and d-hole respectively. The first 
equation from (|^) determine the value of the momentum of the u quark as fc^ = (^ + _E)^ — m^. In the similar fashion 
the second equation from (|^) gives the value of the momentum of the reflected d hole as = (/i — — m^. For z > 
the quasiparticles excitation are described by the wavefunction '0> {t, z) given by the expression ( [10[ ) with B = D = Q 
(no spin flip in the scattering process). Let us note that for < < |A|, the momenta of the quasiparticles has to 
be complex number: 



91,2 = ±VM'-m2-|e|2±2z|e|M (12) 

Particularly simple form of momenta one can find in the massless limit. In that case gi_2 ~ ±/i + i|^|, where the signs 
are chosen as to describe the vanishing of the wavefunction in 2SC phase for z — > cx). Now it is seen that for i? < |A| in 
the 2SC phase the quasiparticle wavefunction is exponentially decaying with the suppression factor |<^| = -^/jAp — E^. 
The continuity conditions matching the wavefunctions are of the form: 

V'< (i, z 0) V> (i, z 0) (13) 
Using this condition one can find the amplitudes of the scattering process: 



E-£,mE f .S\ ^ I 



i.^J|^exp(..) + 0(i 
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In the limit where |A|, i?, m << /i only A and H contribute which means that only hole of the d-green quark can be 
reflected into the QGP phase and the only quasiparticle with momentum qi can propagate in the 2SC phase. It is 
interesting to point out that in the massless limit F — C — exactly. Additional insight we gain from the calculation 
of the probability current. From the equations (|^) it follows that the probability current j = ip^^cttp'fj + ipQ^atpQ is 
conserved. Using ( |lO| , [T^ ) one finds: 

[0 for £; < |A| 



The result ( p^ ) has simple interpretation. If the quark u have energy below the gap it can not excite quasiparticle 
inside the 2SC phase. In that case the matter waves does not propagate through the superconducting medium. 
However if the u quark possesses energy above the gap it excites the quasiparticles with given transition coefRcient. 
This result may be of importance for the physics of the Neutron Stars. 

IV. CONCLUSIONS 

In this paper we describe the phenomena of Andreev reflection at the junction between the cold quark-gluon plasma 
and superconducting 2SC phase. We found that the hole of different flavor and color than the incoming particle is 
reflected toward the QGP phase. Inside the 2SC phase the quasiparicles are excited, whic h for ener gy E of the 
incoming particle below the gap |A|, penetrate the 2SC phase for depth of the order of l/-\/|Ap — whereas for 
the case of i? > |A| they propagate freely inside the superconducting phase. This result is similar to the situation 
one encounters in the condensed matter systems. This can be expected because in the high density QCD the purely 
relativistic phenomena are suppressed by the powers of quark chemical potential. 

From the equation ( p^ one can see that the transport phenomena (like heat transport or density waves) is strongly 
affected by the presence of the interfaces inside the Neutron Stars. The suppression is exponential of the form 
exp (— |A|/r) where T is a temperature of the Neutron Star, usually much smaller than the expected superconducting 
gap. This phenomena certainly require more extensive analysis]^. 

The structure of the Neutron Stars can be complicated, thus more work is needed for better understanding of the 
Andreev reflection phenomena in superconducting QCD. In particular one has to consider other possible supercon- 
ducting phases like CFL or other possible interfaces. Existence of such interfaces certainly influences the dynamics of 
matter flow inside the Neutron Stars. 
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^Let us note that probability current through the interface is not suppressed completely, because in the 2SC phase we have 
unpaired blue quarks. 
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